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ABSTRACT 



. Aims. We have used mid-infrared long-baseline interferometry to resolve the circumstellar material around the Herbig Ae star HR 5999, pro- 
j — ' viding the first direct measurement of its angular size, and to derive constraints on the spatial distribution of the dust. 

Methods. MIDI at the VLTI was used to obtain a set of ten spectrally dispersed (8 - 13yum) interferometric measurements of HR 5999 at 
T^j- different projected baseline lengths and position angles. To derive constraints on the geometrical distribution of the dust, we compared our 
t^- ' interferometric measurements to 2D, frequency-dependent radiation transfer simulations of circumstellar disks and envelopes. 

Results. The derived visibility values between ~0.5 and ~0.9 show that the mid-infrared emission from HR 5999 is clearly resolved. The 
characteristic size of the emission region depends on the projected baseline length and position angle, and it ranges between ~ 5 - 15 mil- 
liarcseconds (Gauss FWHM), corresponding to remarkably small physical sizes of ~ 1-3 AU. For disk models with radial power-law density 
, distributions, the relatively weak but very extended emission from outer disk regions (^ 3 AU) leads to model visibilities that are significantly 
O ! ; lower than the observed visibilities, making these models inconsistent with the MIDI data. Disk models in which the density is truncated at 
' outer radii of ~ 2 - 3 AU, on the other hand, provide good agreement with the data. 
i { Conclusions. A satisfactory fit to the observed MIDI visibilities of HR 5999 is found with a model of a geometrically thin disk that is truncated 
' at 2.6 AU and seen under an inclination angle of 58° (i.e. closer to an edge-on view than to a face-on view). Neither models of a geometrically 
thin disk seen nearly edge-on, nor models of spherical dust shells can achieve agreement between the observed and predicted visibilities. The 
r* ■ reason why the disk is so compact remains unclear; we speculate that it has been truncated by a close binary companion. 

• i-H , 

i Key words, techniques: interferometric - stars: individual: HR 5999 - stars: formation - stars: circumstellar disks 
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1 . Introduction could be equally well explained with such very different mod- 
els as geometrically thin accretion disks and spherical en- 

Herbig AeBe stars (HAEBEs; see Herbig 1960; Herbig 1994; velopes (Hillenbrand et al. 1992; Miroshnichenko et al. 1997). 

The et al. 1994) are pre-main sequence stars of intermediate More complex mode ls, like envelopes with cavities produced 

mass (-2-10 M ) that show broad emission lines, rapid vari- by bipolar mo i eC ular outflows (Hartmann et al. 1993), ge- 

ability, and excess infrared and millimeter-wavelength emis- ome trically thin and optically thick disks surrounded by a 

sion. These properties suggest the presence of substantial spherical enve lope (Natta & Kriigel 1995, Natta et al. 2001; 

amounts of circumstellar dust and gas. The distribution of the Miroshnichenko et al. 1999), flared outer disks, puffed-up in- 

circumstellar material is still a matter of debate. Since until re- ner disk rims (Du ii em0 nd et al. 2001), and disk plus inner halo 

cently the spatial scales of the inner circumstellar environment models ( yinkovic et al. 2006) have also been used to success- 

(a few AU, corresponding to < 0.1") were not accessible to Mly fit me observed SEDs of HAEBEs. The strongest evidence 

optical and infrared imaging observations, conclusions drawn for circumste n ar disks around HAEBEs comes from millime- 

on the spatial distribution of the circumstellar material of the ter i nter ferometry (Mannings & Sargent 1997, 2000): flattened 

HAEBE stars were in most cases based entirely on the model- stmctures a^nd several sources have been resolved on ~ 100 

ing of the spectral energy distribution (SED). However, fits to AU scales at m iiii me ter wavelengths and detailed kinematic 

the observed SEDs are highly ambiguous (e.g. Men'shchikov mode li ng of one source, MWC 480, is consistent with rotation 

& Henning 1994). For example, the SEDs of some HAEBEs m a Keplerian disk 



Send offprint requests to: Thomas Preibisch Speckle interferometry by Leinert et al. (2001) in the NIR 

* Based on observations with the ESO Very Large Telescope reveals numerous examples of resolved halos on scales of 

Interferometer at Paranal Observatory under programme IDs 073C- ~ 0.1"— 1" (corresponding to several hundred AU). Bispectrum 

0248 and 073C-0720. speckle interferometry revealed complex asymmetric struc- 
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tures on scales down to 70 mas in several intermediate- to high- 
mass young stellar objects (e.g. Schertl et al. 2000; Weigelt et 
al. 2002ab, 2006; Preibisch et al. 2001, 2002, 2003; Hofmann 
et al. 2004). Interferometric aperture masking Keck observa- 
tions of the young system LkHa 101 (Tuthill et al. 2001; 2002) 
revealed some details of the inner disk, which is seen nearly 
face-on. To make further progress on this topic, observations 
with even higher spatial resolution, at the level ;$ 10 milli- 
arcseconds (mas), are needed to explore the inner circumstellar 
regions of these objects in detail. Such observations can only 
be performed with long-baseline interferometry. 

The current generation of infrared long-baseline interfer- 
ometers is sensitive to the distribution of dust around the near- 
est young stars on scales on the order of 1 AU, and it provides 
a powerful probe of the models of disks and envelopes of such 
stars. The first systematic study of HAEBEs using the tech- 
nique of long-baseline interferometry in the near-infrared was 
conducted by Millan-Gabet et al. (2001), who could resolve 
the near-infrared emission of 11 of the 15 HAEBEs observed 
at the IOTA interferometer. The derived sizes of the observed 
structures are 0.5 to 6 AU, but no unique geometrical model 
could be determined for any of the sources due to the very 
sparse coverage of the (u, v)-plane. A number of further studies 
have been obtained at different interferometers (e.g., Akeson et 
al. 2000; Wilkin & Akeson 2003; Eisner et al. 2004; Monnier 
et al. 2005). Most of the near-infrared interferometric data are 
consistent with a simple disk model possessing a central dust- 
free cavity, ringed by hot dust, but other geometries can often 
not be excluded, due to the poor sampling of these observations 
of the (u, v)-plane. 

The first long-baseline interferometric observations of 
HAEBEs at mid-infrared wavelengths were recently presented 
by Leinert et al. (2004). Their data, obtained with MIDI at the 
VLTI, trace the 8 - 13 yum emission from hot and warm dust in 
the inner disk regions and can be used to constrain the geomet- 
rical structure of the circumstellar material on angular scales of 
~ 10 mas. They resolved the circumstellar material around all 
of the seven observed, nearby HAEBEs and derived character- 
istic dimensions of the emitting regions at 10 /im ranging from 
1 AU to 10 AU, but, again, they could not determine the geom- 
etry of the dust distributions due to the limited (u, v)-coverage 
of the observations. 

HR 5999 (= HD 144668 = V856 Sco) is one of the best 
studied HAEBEs. It lies in the central part of the Lupus 3 dark 
cloud, which is part of the extended Lupus star-formation com- 
plex. HR 5999 has a spectral type of A7 III-IV and forms a 
~ 45" wide proper-motion binary system with the A1.5 star 
HR 6000. The Hipparcos parallax for HR 5999 corresponds 
to a distance of 208^ pc. The spectral energy distribution of 
HR 5999 shows a clear infrared excess at wavelengths above 
1 iim. IRAS detected HR 5999 at 12, 25, and 60 //m, but pro- 
vided only an upper limit at 100/im. Henning et al. (1994) 
found weak millimeter emission from HR 5999, but the de- 
tection was only at the 1 cr significance level. The far-infrared 
and mm-emission of HR 5999 suggests the total mass of the cir- 
cumstellar material to be on the order 0.006 M (Siebenmorgen 
et al. 2000). HR 5999 belongs to the "group I objects" as de- 



fined by Hillenbrand (1992); the infrared excess of these ob- 
jects is thought to be from emission arising in a flat, optically 
thick accretion disk. The fit of the observed SED with a sim- 
ple analytical model of an optically thick circumstellar disk 
yielded an outer radius of ^ 23 AU and a possible inner hole of 
~ 0. 1 AU. Siebenmorgen et al. (2000) analyzed mid-infrared 
ISO spectra of HR 5999 and found a broad silicate emission 
feature with a peak around 9.6 fim. With a ratio of peak- to 
continuum-flux of about 1 .4, the intensity of this feature is rel- 
atively small. The broad and rather flat shape of the silicate fea- 
ture indicates the presence of relatively large (^ 1 yum) silicate 
grains (cf. van Boekel et al. 2003). 

The stellar parameters of HR 5999 suggest that the object 
is ~ 0.5 Myr old and has a mass of ~ 3-4 M Q (van Boekel 
et al. 2005). HR 5999 will thus arrive at the main sequence as 
a mid B-type star. HR 5999 shows only moderate extinction 
(Ay = 0.49 mag) and is a fast rotator (vsin/ = 204 km/sec; 
Royer et al. 2002). The quasi-periodic and irregular photomet- 
ric and spectroscopic variability of HR 5999 in the optical (see 
Perez et al. 1992) and in the UV range (see Perez et al. 1993) 
has been interpreted as due to instabilities in an optically thick 
circumstellar accretion disk. The detailed analysis of the UV 
emission lines in the IUE spectra of HR 5999 by Perez et 
al. (1993) revealed gas accreting onto the star with velocities 
as high as 300 km/sec (see also Blondel et al. 1993). Perez et 
al. (1993) also concluded that the observed double-peaked Ha 
and Mg II emission profiles strongly suggest that the circum- 
stellar disk around HR 5999 is seen nearly edge-on. The aver- 
age mass accretion rate was estimated to be ~ 7 x 10~ 7 M /yr, 
but the observed variability of the UV emission clearly suggests 
that the accretion process happens in a non-steady fashion. 

To summarize, the existing data clearly show that HR 5999 
is surrounded by circumstellar material, probably in the form 
of an accretion disk. The high vsin / and the double-peaked 
structure of the UV emission lines suggest that the system is 
seen under a relatively high inclination angle (i.e. close to edge 
on). This makes it an interesting target for a detailed interfero- 
metric study of the geometry of the circumstellar dust. Grady 
et al. (2005) presented HST STIS white-light coronographic 
imaging results for HR 5999, found no indication of nebulosity 
at r > 70 AU, and concluded that the disk around HR 5999 
must be geometrically small. 

HR 5999 has a close visual companion, Rossiter 3930 
(see Stecklum et al. 1995, and references therein), with AV ~ 
4.6 and AK ~ 3.6, at a projected separation of 1.4", cor- 
responding to a projected distance of about 300 AU. Since 
Rossiter 3930 shows strong Ha emission, it is very likely a 
T Tauri star. If Rossiter 3930 is actually physically related to 
HR 5999 and not a chance projection, it may be responsible for 
the rather low (sub-)mm emission and corresponding low mass 
of the circumstellar material of HR 5999 due to the mecha- 
nism of dynamic disk clearing by gravitational binary interac- 
tion (e.g. Arymowicz & Lubow 1994). 

2. Observations and data analysis 

The observations discussed in this paper were obtained with 
the MIDI instrument at the VLTI during several observing runs 
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Table 1. Log of MIDI observations of HR 5999. 



UT date & time 


B 


PA 


Obs. 






[m] 


[deg] 


name 


2004-04-10 


02:44:34 


46.34 


173.37 


46m/173° 


2004-04-10 


05:23:56 


46.60 


21.57 


46m/21° 


2004-04-11 


05:16:58 


46.59 


20.74 


46m/20° 


2004-04-11 


09:09:04 


42.42 


52.12 


42m/52° 


2004-04-12 


09:46:20 


39.76 


56.36 


39m/56° 


2004-06-28 


00:15:43 


102.38 


15.35 


102m/15° 


2004-06-28 


02:06:06 


100.76 


30.90 


100m/30° 


2004-06-28 


04:26:08 


90.89 


46.10 


90m/46° 


2004-06-28 


05:22:28 


83.22 


50.57 


83m/50° 


2004-09-30 


00:02:18 


63.50 


123.74 


63m/123° 



Table 2. List of calibrators used for the calibration of the 
HR 5999 data together with uniform-disk diameters (c/ud) and 
the date, as well as the time stamp (? bs) of MIDI observations. 



Calibrator 
HD number 




(mas) 


night 
(2004) 


t obs (UTC) 


81797 


9.54 


±0.75 


2004-04-10 
2004-04-11 


00:41:19 
03:56:19 


107446 


4.54 


±0.23 


2004-04-10 
2004-04-11 
2004-04-12 


02:13:31 
23:55:03 
01:24:27 


109379 


3.30 


±0.17 


2004-06-27 


23:27:48 


129456 


3.39 


±0.09 


2004-04-10 
2004-04-12 
2004-04-12 
2004-06-28 


05:46:05 
08:37:42 
10:09:08 
00:43:13 


134505 


2.55 


±0.13 


2004-06-28 


01:43:12 


139127 


3.30 


±0.22 


2004-04-11 
2004-04-11 
2004-04-12 
2004-06-28 
2004-06-28 


05:35:48 
09:36:42 
07:09:25 
02:48:50 
05:43:29 


152885 


2.88 


±0.09 


2004-04-10 
2004-04-12 
2004-04-12 


07:13:48 
05:57:26 
09:24:41 


161892 


3.94 


±0.21 


2004-04-12 
2004-04-12 


04:42:03 
05:02:42 


168723 


2.88 


±0.13 


2004-04-10 
2004-04-10 
2004-09-30 


09:29:56 
10:12:13 
01:22:01 


169916 


3.39 


±0.09 


2004-04-12 


07:57:12 


176411 


2.24 


±0.13 


2004-04-10 
2004-04-10 


08:08:34 
08:45:41 


188512 


2.07 


±0.10 


2004-09-30 


03:13:45 



between April and September 2004. In total, ten MIDI mea- 
surements of HR 5999 were performed: six visibility measure- 
ments for observing program 073C-0248 (PI: Preibisch) and 
four visibility measurements for program 073C-0720 (PI: van 
Boekel). The combined data set is analyzed in this paper. The 
observation dates, the projected baseline lengths, the position 
angles of the projected baseline on the sky, and the abbreviated 



Fig. 1. Points in the (u, v) plane sampled by our MIDI observa- 
tions of HR 5999. 



data set names used in the following text are summarized in 
Tabled Our observations cover a range of projected baseline 
lengths from 39 m to 102 m and position angles ranging from 
15° to 173°, and, therefore, provide a relatively good coverage 
of the (m, v) plane. 

The fringes were scanned by applying rapid internal optical 
path difference (OPD) changes of a few wavelengths around 
the point of zero OPD. A prism with a spectral resolution of 
A/AA ^ 30 was used to obtain spectrally dispersed fringes 
in the range of 8 to 13 /im. Immediately after a recording 
of interferometric data, photometric calibration data were ob- 
tained by blocking each of the two beams alternately. Chopping 
with a frequency of 2 Hz was applied to obtain the data 
for the sky background subtraction. A more detailed descrip- 
tion of the observing procedure is given in Przygodda et al. 
d2003l . |Leinert et al. (2004)| Chesneau et al. (20051 I2005b> . 
and Ohnaka et al. (I2005i . 

Calibration stars with known angular diameter were ob- 
served for the calibration of the raw visibilities. Table Ogives 
an overview of these calibration stars, their uniform-disk di- 
ameters as given in the CalVin list 1 , and the MIDI observa- 
tions that we used for the calibration of the raw visibilities of 
HR 5999. 

We utilized two different MIDI data reduction packages: 
MIA developed at the Max-Planck-Institute for Astronomy and 
EWS developed at the Leiden Observatory. While the MIA 
package is based on the power spectrum analysis and measures 
the total power of the observed fringes, the EWS package first 
corrects for OPDs (instrumental, as well as atmospheric de- 
lays) in each scan, and then the interferograms are averaged. A 
more detailed description of the MIA and EWS packages can 
be found in Leinert et al. (2004) and Jaffe (2004l, respectively. 
In our analysis of the HR 5999 data, we found good agreement 
between the results of both packages within the error bars. In 



available from ESO at http://www.eso.org/observing/etc/ 
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the following discussion we always refer only to the results ob- 
tained with the MIA software. 

The errors of the calibrated visibility were estimated from 
the 1 <x scatter of the visibilities calibrated with all individual 
data sets of a night. In addition, to account for systematic error 
contributions (e.g. due to imperfect beam overlap), we assumed 
a minimum visibility error of 15%. One of the ten visibility 
measurements, the 46m/21° data set, was affected by strong 
variations in the transfer function around the time of the mea- 
surement. This particular data set could not be calibrated well; 
we decided to drop it from the following analysis. 

3. Observed visibilities 

As the first step of the analysis, we consider the dependence of 
the observed visibilities on the wavelength, the projected base- 
line length, and the baseline position angle. Figure[5]shows the 
wavelength-dependence of the visibilities of HR 5999. All vis- 
ibility curves are rather flat, i.e. the observed visibilities show 
no strong dependence on the wavelength. 

Figure|3]shows the observed 10.5 pm visibilities as a func- 
tion of projected baseline length. One can see that the data are 
not consistent with simple uniform-disk or Gaussian brightness 
distribution models. 

3.1. Size of the emission region 

In order to derive an initial rough estimate of the size of the re- 
solved emission region, we fitted each observed visibility value 
with two simple models: first, a uniform-disk model, and sec- 
ond, a Gaussian brightness distribution. The resulting diame- 
ters range between 10 and 25 milliarcseconds (corresponding 
to ~ 2 - 5 AU at the distance of 208 AU) for the uniform 
disk model and 5-15 mas (~ 1 - 3 AU) for the Gauss model 
(FWHM). The Gauss-fit diameters as a function of wavelength 
for the different MIDI observations are shown in Fig. [4] For all 
data sets, the diameters increase with increasing wavelength in 
the MIDI band. This effect can be understood as a consequence 
of the fact that the emission at longer wavelengths comes from 
cooler material, which is located at larger distances from the 
central star than the warmer material radiating at shorter wave- 
lengths. Typically, the diameters at 13 fim are ~ 50% - 100% 
larger than at 8 /im. 

Another notable effect is the trend that the longest baselines 
give the smallest diameters. Although one must not ignore that 
the different observations were obtained at different position 
angles (see next paragraph), this suggests a complex geometry 
for the emission region, as discussed in the next section. 
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Fig. 2. MIDI visibilities of HR 5999 plotted versus wave- 
length. The upper panel shows the data for the observations at 
baseline-length/ position angles of 102m/15° (dots), 83m/50° 
(crosses), 63m/123° (boxes), and 46m/20° (diamonds). The 
lower panel shows the data for the observations at 100m/30° 
(dots), 90m/46° (crosses), 46m/173° (boxes), 42m/52° (dia- 
monds), and 39m/56° (asterisks). The thick dash-dotted line 
and the corresponding scale on the right y-axis in the upper 
panel show the flux-calibrated spectrum obtained with MIDI 
during the 39m/56° observation. 



3.2. Shape of the emission region 

In Fig.[5]we plot the Gauss-fit radii (at wavelength 10.5 /im) as 
a function of the position angle of the observation. The plot 
suggests that the structure of the emission region is not de- 
scribed well by a spherically symmetric distribution, i.e. seems 
to be asymmetric. An inclined ellipse with axis ratio 1:1.56 ap- 
parently yields better agreement to the data than a circle, but 
does not provide a particularly good fit either. The data suggest 



that the emission has a more complicated structure, probably 
with substructure at different size scales. A more detailed anal- 
ysis of size and shape of the emission region is described below 
in Section|4] 
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Fig. 4. Gauss-fit FWHM diameters as a function of wavelength derived for the MIDI observations of HR 5999. The individual 
observations are identified by their values for projected baseline length and position angle. 
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Fig. 3. Visibility at 10.5 fim as a function of projected baseline 
length. The data points are annotated with their respective po- 
sition angles. For comparison, we also plotted the theoretical 
visibility curves for uniform disk models with diameters of 10, 
15, 20, and 25 mas (solid lines) and for Gauss models with 
FWHM = 5, 10, and 15 mas (dashed lines). 



3.3. Comparison to MIDI data of other HAEBEs 

The visibilities measured for HR 5999 differ from those ob- 
tained by Leinert et al. (2004) for a sample of seven HAEBEs in 
two ways. First, all MIDI measurements of HR 5999 give visi- 



Fig.5. Gauss-fit FWHM diameters at 10.5 yum for the differ- 
ent observations as a function of baseline position angle. The 
dotted circle has a radius of 3.7 mas. The dashed ellipse has a 
minor axis of 4.8 mas and a major axis of 7.5 mas. 

bilities > 0.5, whereas six of the seven HAEBEs in the Leinert 
et al. (2004) sample show much lower visibility values, < 0.4 
in the 8yum to 13/Lmi wavelength range. Only 51 Oph shows 
visibility values similar to what we find for HR 5999. This im- 
plies that the angular extent of the mid-infrared emission from 
HR 5999 is smaller than that of most HAEBEs in the Leinert 
sample. 
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The second difference is the wavelength dependence of the 
visibilities. While the visibility for HR 5999 is essentially flat 
over the 8 yum to 13 yum wavelength range, most of the HAEBE 
stars in the Leinert sample show a rather strong drop in the vis- 
ibility between 8^m and 10 //m, in some cases followed by a 
slight increase towards longer wavelengths. The latter wave- 
length dependence of the visibility is expected for objects with 
significant 10 fim silicate emission features (a strong emission 
feature in the spectrum produces a minimum in the visibility 
due to optical depth effects, because the apparently larger ob- 
ject size in the center of the emission line yields lower visibili- 
ties). The flat visibility curves of HR 5999 are probably related 
to its very weak silicate emission feature, which does not con- 
tribute enough to produce a significant dependence of the visi- 
bility (apparent object size) on wavelength over this feature. 

4. Constraints on the source shape from 2D 
radiative transfer simulations 

The true morphology of the mid-infrared emission is most 
likely more complex than simple geometrical models for the 
brightness distribution, like uniform-disk or Gauss models. For 
example, the existence of a dust-free hole inwards of the dust 
sublimation radius already implies that the brightness distri- 
bution cannot have a steadily decreasing profile. In order to 
see how much information on the spatial distribution of cir- 
cumstellar material can be deduced from our MIDI data, we 
performed 2D radiation transfer calculations. We believe that 
this approach provides more realistic constraints than the use of 
simple analytical models. Our radiative transfer modeling also 
includes a self-consistent determination of the dust sublimation 
radius, i.e. the inner disk edge, which is a crucial ingredient in 
any realistic modeling attempt. Note that the aim of the mod- 
eling in this paper is an exploration of different geometries for 
the dust distribution and not a detailed disk modeling (which is 
not warranted considering the rather large uncertainties of the 
measured visibilities). 

4.1. Modeling procedure 

For our simulations we used the 2D radiation transfer code 
described in detail by Sonnhalter et al. (1995). In this code, 
the distribution of dust temperatures and radiation intensi- 
ties in an axially symmetric dusty circumstellar environment 
around a central radiation source is calculated within the frame- 
work of the flux-limited diffusion approximation (Levermore & 
Pomraning 1981). The equations are discretized on a quadratic, 
equally spaced grid. As the density distribution is assumed to 
be axially symmetric with respect to the z-axis and mirror sym- 
metric with respect to the plane z — 0, cylindrical coordinates 
(r, z, 0) are used, because the symmetry reduces the relevant 
information to the first quadrant of the r - z plane. To im- 
prove resolution and convergence, a system of 5 nested grids 
with decreasing grid spacing is used. The differential equations 
for the radiation field are iterated together with the tempera- 
ture equations until a self-consistent equilibrium configuration 
is reached. At the end of the iteration procedure, the total radia- 
tion flux (integrated over all frequencies) through the outermost 



grid cells was always within 1% of the stellar luminosity for 
all computed models. We note that any non-spherical density 
distribution results in non-spherically symmetric temperature 
distributions. For example, in models of a circumstellar disk, 
the disk surface is generally found to be warmer than the disk 
midplane at the same radial distance from the star. 

After determination of the dust temperatures for a multi- 
component dust model, a ray-tracing procedure is used to cal- 
culate intensity maps for the appearance of the central object 
and its circumstellar environment at varying inclinations for 
selected frequencies. The simulated images are then used to 
compute visibility maps in the (u, v) plane, which finally can be 
compared to the MIDI measurements. We used 64 wavelength 
points between 0.1 fim and 5 mm for our modeling of HR 5999. 
The size of the numerical grid was set to 50 x 50 AU, which re- 
sulted in a spatial resolution (i.e. the grid point spacing for the 
innermost grid) of 0.016 AU. As a final check on the numerical 
accuracy, we used images computed with the ray-tracing pro- 
cedure at a large number of different inclination angles to deter- 
mine the total flux (integrated over all wavelengths) as a func- 
tion of the inclination angle. Integration of these fluxes over all 
inclination angles then gives the total model luminosity, which 
finally can be compared to the input luminosity. For the geo- 
metrical models described below, we find typical deviations of 
^ 1 % in the case of spherical and thick disk models, increas- 
ing up to ~ 8% for the thinnest disk models (when the vertical 
thickness of the inner disk gets similar to the numerical grid 
resolution). 

Note that our radiative transfer code not only accounts for 
extinction by and emission from dust grains, but also treats 
scattering by dust grains; isotropic scattering is implicitly as- 
sumed. Although the scattering efficiency of dust grains is quite 
low at mid-infrared wavelengths, scattering of optical/UV pho- 
tons in the region above and below the disk plane can lead to 
significantly increased heating of the outer disk regions; this 
can result in more extended thermal emission at mid-infrared 
wavelengths. If scattering is neglected (as is sometimes done 
in other radiation transfer simulations), the extent of the emis- 
sion can be significantly underestimated. 

Our radiative transfer modeling assumes the disk is passive 
(i.e. the temperature at each point is determined by the local 
radiation field) and does not take into account possible viscous 
disk heating by accretion. This is appropriate because the stel- 
lar parameters and the estimated accretion rate (see § 1) suggest 
that accretion does not contribute more than ~ 10% to the total 
luminosity. 

A self-consistent radiation transfer model should, of 
course, reproduce not only the observed visibilities, but also the 
broad-band spectral energy distribution (SED). We therefore 
constructed the SED of HR 5999 from the UV to mm wave- 
lengths from literature data (Hillenbrand et al. 1992; Henning 
et al. 1994; Acke & van den Ancker 2004). A complication 
arises from the presence of the companion Rossiter 3930, 
which has an angular separation of only 1.4" from HR 5999. 
The fluxes for the SED points, which were obtained from 
seeing-limited observations, therefore represent the sum of 
the fluxes from HR 5999 and Rossiter 3930. The adap- 
tive optics study by Stecklum et al. (1995) showed that the 
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flux of Rossiter 3930 is considerably smaller than that of 
HR 5999 at optical and near-infrared wavelengths; they de- 
termined A(V,J,H,K) ~ (4.6, 3.0, 3.1, 3.6) mag. At longer 
wavelengths, however, Rossiter 3930 may well contribute a 
larger fraction to the total flux, and it cannot even be excluded 
that it may dominate the total flux, as is the case for the sev- 
eral well-known "infrared companions" in binary T Tauri star 
systems (e.g. Chelli et al. 1988). Although the latter possibil- 
ity does not appear to be very likely in the case of HR 5999 
(since there is no strong increase in the flux contribution from 
Rossiter 3930 from the optical to the near-infrared range), all 
SED points longward of 2.2 /mi can only be considered as up- 
per limits to the flux of HR 5999. In our radiation transfer mod- 
eling we therefore cannot require that the model SED matches 
the observed SED longward of 2.2 fim, but only that it does not 
exceed the observed SED. We note that the companion does 
not affect our MIDI data of HR 5999 since the interferometric 
field-of-view is smaller than the separation 2 . 

As input for the radiation transfer calculations one has to 
specify the luminosity and effective temperature of the central 
radiation source (we used L^ \ = 87 L and r e ff = 7925 K fol- 
lowing van Boekel et al. 2005), the optical properties of the 
dust, and the density structure of the circumstellar material. 
The dust composition around HR 5999 was recently investi- 
gated by van Boekel et al. (2005), who found that the 8-13 pm 
spectrum is dominated by large silicate (olivine) grains. Guided 
by their results, we used a dust model here that consists of sili- 
cate and carbon grains with sizes of 1 pm. The sublimation tem- 
perature of both grain species is taken to be 1500 K. In the radi- 
ation transfer code, the destruction of grains through sublima- 
tion is simulated by setting the number density of the grains to 
zero, if the temperature of that particular grain species, which 
is calculated simultaneously with the radiation field, exceeds 
the sublimation temperature. 

Model images were computed with the radiation transfer 
code for those 15 wavelengths between 8yum and 13 fim for 
which the visibilities were determined from the MIDI data. 
We then calculated visibility maps from the model images and 
compared them to the measured visibilities for HR 5999. Since 
the orientation of the disk (i.e. the angle between the disk axis 
and the north-south direction on the sky, measured counter- 
clock wise and denoted here as the position angle <p) is un- 
known, an x 1 fitting procedure was used to determine the po- 
sition angle that results in the best agreement between mea- 
sured and simulated visibilities. For this we used all 9 visibility 
curves, each consisting of 15 wavelengths points, i.e. 135 indi- 
vidual data points. 



4.2. Models for the density distribution 

We consider here two different classes of models for the den- 
sity distribution of the circumstellar material, disk models, and 
spherical envelope models. 



4.2.1 . Extended disk (Model D) 

We started with models of a quasi-Keplerian disk (i.e. the disk 
structure is determined by the gravitationally dominating cen- 
tral star), with the following density distribution: 



Pd 



/>o| — | exp 



n I z 



(1) 



with h z = \ — \ x h- r . 
Vol 

We used p = 15/8 for the radial density power-law index, 
and for the radial dependence of the vertical disk scale height 
h z we assumed q = 9/8, i.e. a slightly flaring disk shape. The 
parameter h fixes the relative geometrical thickness of the disk. 

In our modeling we considered different values for the rela- 
tive disk thickness h (from very thin disks with h = 0.05 to very 
thick disks with h ~ 0.8), different inclination angles i (from 
i = 0° for a face-on view and i = 90° for an edge-on view), and 
also tried different values for the radial density power-law index 
p. For each combination of model parameters, the density po is 
restricted by the requirement that the model SED should agree 
with the observed SED points for wavelengths A < 2.2 pm and 
must not exceed the observed fluxes for longer wavelengths 
(where the observed SED points contain an unknown flux con- 
tribution from the binary companion). 

The resulting brightness distribution in the simulated im- 
ages is dominated by the central star and the inner disk region 
of a few AU (see Fig.|6j. The peak brightness is found at the 
inner rim of the dusty disk at the dust sublimation radius, and 
it drops smoothly outwards. Despite the quite compact appear- 
ance of the simulated images, the weak emission from the outer 
parts of the disk contributes a significant fraction of the total 
flux at mid-infrared wavelengths. This extended emission leads 
to a large effective size of the emission and yields model visi- 
bilities that are considerably lower than some of the observed 
visibilities. The characteristic size of the model brightness dis- 
tributions is up to twice as large as the sizes derived from the 
observed visibilities. Even if we change the power-law expo- 
nent for the radial density law (to values as steep as p = 4) or 
the thickness of the disk, the extent of the emission in the model 
images is still considerably too large and inconsistent with the 
observed visibilities. 

4.2.2. Truncated disk (Model TD) 

One possibility for decreasing the amount of extended emis- 
sion in the outer parts of the images is to limit the circumstellar 
mass distribution to a certain outer radius. For this, we used the 
disk model describe above (Equ.Q and truncated the density 
distribution at the outer radius r,. This truncation is performed 
by multiplying the density distribution with a Fermi-type func- 
tion: 



Ptd = Pd x 1 1 + exp I 



(2) 



2 To be sure, we inspected the acquisition images and found no in- 
dication of another source in the MIDI field-of-view. 



where po is the density distribution given in Equ. Q, R = 
y r 2 + z 2 j s me rac [ial distance from the origin, while a deter- 
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Fig. 6. Radiation transfer modeling of HR 5999. The upper row shows the model images at lOjjm for the disk model (D, left), 
the truncated disk model (TD, middle), and the spherical shell model (S, right). Each image shows a 12.5 x 12.5 AU region with 
a negative logarithmic intensity scale, seen under an inclination angle of i = 59° and at position angle = 24°. The mosaic of 
visibility versus wavelength plots compares the observed visibilities of HR 5999 (solid dots with error bars) for each individual 
MIDI observation to the model visibilities; the truncated disk model (TD) is displayed by the solid line, the disk model without 
truncation (D) by the dashed line, and the spherical shell model (S) by the dotted line. 



mines the width of the transition zone of the Fermi-type func- 
tion; we used a = R,/1Q. 

With these models we were able to find acceptable fits to 
the observed SED up to A ~ 10 fim (see Fig. [7J and the ob- 
served visibilities for truncation radii in the range ~ 2 - 3 AU, 
inclination angles of i ~ 60°, and position angles of <p ~ 25°. 
The vertical thickness of the disk (governed by the parameter h) 
is nearly unconstrained by the data (thin and moderately thick 
disks seen under intermediate inclination angles "look" very 
similar). 



In Fig. [6] we show a model with parameters h = 0.056, 
R, = 2.6 AU, i = 59°, and <j) = 24°, which yields an acceptable 
(albeit not perfect) fit to the data with;^ = 1.07. The V-band 
optical depth of this model along the disk plane is r = 3.5. 
Although models with larger vertical disk scale height give 
marginally better fits to the visibilities (for example, x 2 r - 0.98 
for a disk with h = 0.4), we prefer this specific thin-disk model 
because the vertical scale height of the dusty disk agrees with 
the formal condition of vertical hydrostatic equilibrium (see, 
e.g., Equ. 7. in Chiang & Goldreich 1997). 



Preibisch et al.: VLTI MIDI observations of HR 5999 



9 




1.0 10.0 100.0 

Wavelength [/^m] 



Fig. 7. Observed spectral energy distribution (SED) of HR 599 
(boxes; compiled from the references listed in the text) com- 
pared the predicted SED for the truncated disk model described 
in the text. While the model reproduces the SED reasonably 
well up to A ~ 10 /mi, it does not account for the bulk of the 
cool dust responsible for the emission at > 13yum. This cool 
dust may be associated to the companion star Rossiter 3930 
or perhaps widely distributed in the HR 5999/R 3930 system. 
(Note that dust at distances > 100 AU from HR 5999 would 
be cool enough to produce excess emission only longward of 
10 pm. The beamsize of the far-infrared SED points corre- 
sponds to sizes of several 1000 AU, and the data do therefore 
not allow to constrain the location of the cool dust.) 



For comparison, we note that a model for the same disk 
parameters but without the truncation factor (model D in Fig. [(J 
yields^ = 18.3. 

4.2.3. Spherical shell (Model S) 

As already shown in 33.21 the observed visibilities suggest 
that the brightness distribution is not spherically symmetric but 
rather somewhat elongated. In order to assess quantitatively 
how much the truncated disk model provides a better descrip- 
tion to the data than spherical shell models, we also considered 
models with a spherically symmetric power-law density distri- 
bution, in which the density distributions were again truncated 
at radius R,: 



Ps = Pi) 



1 + exp 



R-R, 



(3) 



We found that spherical shell models always provide poorer 
fits to the data than the truncated disk models. A spherical 
model (model S in Fig. [6} with the same radial density power- 
law and truncation radius as in the truncated disk model de- 
scribed above yields^ = 2.91. 



4.3. Conclusions from the radiative transfer modeling 

We find that dust distributions in which the radial density fol- 
lows a power law (with exponents in the range p ~ 1-4) 
produce brightness distributions that are much too extended to 
be consistent with our set of MIDI data. Disk models in which 
the density is truncated at outer radii of ~ 2-3 AU, on the 
other hand, provide reasonably good fits to the MIDI data. 

Given the large number of free parameters in the radiative 
transfer models and the relatively large visibility errors, which 
do not allow us to constrain model parameters very well, it was 
not our aim to scan the full multi-dimensional parameter space 
in order to find the overall "best fit" model. We instead pro- 
ceed the other way round and summarize here those model cat- 
egories that can be ruled out by our data and modeling. The 
following cases are inconsistent with the MIDI data: 

- Spherical shells - The observed visibilities suggest that 
the brightness distribution is not circular symmetric but 
rather somewhat elongated. This is confirmed by our radia- 
tive transfer modeling, which shows that spherically sym- 
metric shell models cannot produce good fits to the ob- 
served data. 

- Disks seen at low inclination angles (i < 30°, i.e. nearly 
"face-on") - Such models produce a nearly circular sym- 
metric brightness distribution and thus do not fit the data 
well. 

- Very thick disks (/; > 0.7) - Such models also pro- 
duce only very slightly elongated, nearly circular symmet- 
ric brightness distributions. 

- Geometrically thin disks (h < 0.2) seen close to edge-on 

(z > 75°) - Such models predict strongly elongated bright- 
ness distributions, which are inconsistent with the moderate 
amount of asymmetry suggested by the MIDI data. 

Finally, we note that there may be other ways to explain the 
compact brightness distribution. One possibility for decreasing 
the amount of extended emission in the outer parts of the im- 
ages may be to introduce a dense puffed-up wall at the inner 
edge of the dust distribution that casts a shadow on the outer 
parts of the disk. Such a shadow leads to a much steeper tem- 
perature profile in the disk regions behind the rim, and there- 
fore the region containing warm enough dust that radiates sig- 
nificantly in the mid-infrared is smaller than without an in- 
ner rim. Models of this kind have, for example, been used by 
Dullemond et al. (2001) to fit the observed SEDs of HAEBEs. 
The resulting brightness distribution will depend on the de- 
tailed shape of the inner rim, which is not well known (see, e.g., 
discussion in Isella & Natta 2005). Since such models intro- 
duce several additional free parameters for the radiative trans- 
fer modeling, they are beyond the scope of the present study 
and will be the topic of a forthcoming paper. 

5. Summary and conclusions 

Our MIDI interferometric study has, for the first time, resolved 
the mid-infrared emission around the Herbig Ae star HR 5999. 
The most remarkable result is the compactness of the 10 pm 
emission, which has a characteristic radius of only ~ 2 - 3 AU. 
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The data also suggest that the emitting region is elongated. Our 
radiative transfer modeling shows that a moderately thick disk 
with outer edge at ~ 2.6 AU, seen under an intermediate in- 
clination angle, provides a reasonably good description of the 
data. We also find that disk models provide better fits to the 
data than spherical shell models. 

The derived outer disk radius of ~ 2.6 AU is consider- 
ably smaller than canonical assumptions of disk sizes around 
young stellar objects. This leaves us with the question of why 
the disk around HR 5999 is so compact. One possible reason 
may be that the disk was truncated by the gravitational effect 
of a companion to HR 5999. Dynamical disk clearing is ex- 
pected to truncate the disks around the stars in a binary sys- 
tem to roughly 1/3 of the orbital separation of the components 
(see, e.g., the numerical simulations by Pichardo et al. 2005). 
The presence of the companion Rossiter 3930 is interesting in 
this context, but the quite large projected separation between 
HR 5999 and Rossiter 3930 of ~ 300 AU would require an 
extremely elliptical orbit in order to explain a disk truncation 
radius of only ~ 3 AU. A more likely possibility would be the 
presence of another, yet undetected, companion at close sepa- 
ration (~ 10 AU) from HR 5999. The MIDI data provide no 
hint of the existence of a close companion, but we note that it 
would essentially be undetectable in the MIDI data if its mid- 
infrared emission 3 is less than a few percent of that from the 
material around HR 5999. 

Acknowledgements. We would like to thank the referee for detailed 
comments that helped to improve this paper. 
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